Frontiers in Molecular Biosciences

In our past study on human cilia from primary ciliary dyskinesia (PCD) patients using cryo-electron tomography and mass spectrometry, defects in the DNAH5 gene were proven to cause the loss of the entire outer dynein arm along with other proteins. This phenomenon was not observed in the unicellular green alga Chlamydomonas. In this study, we examined the loss of ciliary components caused by defects in various genes, based on proteomic and structural analyses of motile cilia from multiple PCD patients. The absence of single outer arm dynein genes causes a loss or significant decrease in other inner and outer arm dyneins. Interestingly, some patients with outer dynein defects show decreased intraflagellar transport (IFT) proteins, suggesting an influence of cargo components or assembly on transport. These phenomena were not observed in Chlamydomonas, suggesting a more complex mechanism of ciliogenesis in humans despite the high similarity in the final 3D architecture of the 9+2 axoneme. Defects in certain central pair proteins unexpectedly reduced components on the doublets as well. Our results demonstrate protein-protein interactions during stages of human motile ciliogenesis that are more complex than those in unicellular organisms.

Bordetella pertussis, the causative agent of whooping cough, is a reemerging public health threat. While the Tripartite Tricarboxylate Transporter (TTT) system BctCBA was previously implicated solely in citrate uptake, we demonstrate that the solute-binding protein BctC specifically binds citrate chelated with Zn{superscript 2} and Ni{superscript 2}. To elucidate the molecular mechanism of this interaction, we determined the crystal structures of BctC in three states: apo, open, and closed (citrate-zinc-bound), defining the structural determinants for metal-citrate recognition. Comparative analyses suggest that citrate-mediated divalent cation binding is a widespread feature among bacterial TTT homologs. Finally, in silico modeling of the full BctCBA complex predicts an elevator-type transport mechanism. Together, these findings redefine the functional scope of BctCBA, revealing a sophisticated strategy by which B. pertussis exploits organic chelators to acquire essential trace metals during infection.

The limited understanding of transporter substrate spectra constrains our ability to interpret cell and membrane function, highlighting the need for methods that enable transporter deorphanization and characterization of promiscuous transport activities. Here, we present a Xenopus oocyte-based platform for unbiased transporter substrate discovery. Oocytes expressing heterologous solute carrier proteins (SLC) were incubated in human blood serum, a chemically complex metabolite library containing thousands of endogenous metabolites and xenobiotics, followed by paired untargeted LC-MS/MS profiling of intracellular extracts and surrounding medium to capture metabolite exchange events. Across the five human SLC transporters, viz. SLC10A2, SLC10A6, SLC13A2, SLC16A10, and SLC46A1, metabolite exchange signatures were detected, and automated feature annotation was refined by manual chromatographic peak inspection. The workflow recovered known substrates of SLC10A2 and SLC16A10 and identified additional transported metabolites with MS/MS confirmation. This method provides a scalable framework for transporter substrate profiling and prioritization of candidates for targeted validation.

Amylose contributes to starch crystallinity, but the stability of packed amylose double helices in water at elevated temperature remains insufficiently characterized. Here, we used molecular dynamics simulations to test whether chain length affects the short-timescale stability of A-type amylose oligomers in water. Six systems differing in chain length (6, 12, or 24 glucose units per chain) and oligomer size (isolated double strand or dodecamer of six double strands) were simulated, and five independent 1-s production runs were analyzed for each simulated condition. Oligomers with six glucose units showed structural collapse accompanied by increased water penetration. By contrast, dodecamers with 12 or 24 glucose units largely retained packed double-helical organization over the simulated timescale, although fraying was observed at their ends. These results indicate that chain length and lateral packing strongly affect the early structural response of amylose-like crystalline segments in hot water. The present simulations do not establish the ultimate fate of longer oligomers at longer timescales, but they identify a relative stability difference that is relevant to molecular interpretations of hydration-driven disordering in starch.

Nucleotide modifications of the tRNA anticodon can affect protein translation fidelity and speed. Chemical modifications of the anticodon nucleotide 34 are regulated under cellular stress and associated with several translational defects and pathologies. Here, we investigate how these modifications influence A-site codon recognition interactions and their coupling to the CAR site that lies adjacent to nucleotide 34 in the ribosome. The conserved three-residue CAR interface hydrogen bonds in a sequence-dependent manner to the mRNA +1 codon 3-adjacent to the A-site codon and is implicated in tuning translational speed. The C of CAR is pi-stacked with the nucleotide 34 of the A site tRNA anticodon. The A site and the CAR site influence each others hydrogen bonding and stacking interactions, and these codon-adjacency effects potentially provide a layer of regulation affecting translational fidelity and kinetics. Through molecular dynamics simulations of a subsystem of a translocating ribosome IRES-model, we observed that nucleotide 34 modifications affect the hydrogen bonding and stacking interactions at the A site and CAR site as well as CARs influence on the A site interactions. Integrating these results with gene sequence and ribosome profiling analyses, we propose that nucleotide 34 modifications help modulate CARs sequence-dependent tuning of translation in response to cellular stress. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/740946v1_ufig1.gif" ALT="Figure 1"> View larger version (20K): org.highwire.dtl.DTLVardef@fb5a2eorg.highwire.dtl.DTLVardef@1dda97eorg.highwire.dtl.DTLVardef@1a29441org.highwire.dtl.DTLVardef@b353f0_HPS_FORMAT_FIGEXP M_FIG C_FIG

Engulfment and Cell Motility protein 1 (ELMO1) regulates cell migration, phagocytosis, and cytoskeletal remodeling, positioning it as a compelling therapeutic target across kidney diseases, oncology, enteric infections and inflammation. Despite this potential, no approved therapeutics or clinically validated small-molecule modulators of ELMO1 currently exist. ELMO1 functions by forming a complex with DOCK180 (or DOCK2) to activate the small GTPase Rac1, and the recent structural resolution of the ELMO1/DOCK2 complex now provides an opportunity to target this protein-protein interface directly. Here, we present the first investigation into the druggability of the ELMO1/DOCK2 complex and report the initial virtual screening to identify small-molecule inhibitors of this interaction. Molecular dynamics (MD) and free energy level (FEL) studies were carried out to validate the potential of the predicted hits. This work establishes a computational framework for the development of the first generation of ELMO1-targeted therapeutics. In addition to demonstrating the drugability of ELMO1, this work introduces two open-source Python tools for the rapid analysis and visualization of protein-protein interaction and ligand-protein MD trajectories from DESMOND output files. These tools are designed to be broadly accessible, offering practical utility to the wider DESMOND user community.

Eukaryotic translation initiation is tightly regulated by interactions among translation initiation factors (eIFs) that ensure accurate start codon selection. The translation regulator, eIF5 mimic protein 1 (5MP1) contributes to this process by competing with eIF5 for binding to eIF2, thereby increasing the stringency of translation initiation. Despite its important regulatory role and emerging involvement in tumorigenesis, structural information on human 5MP1 remains limited. Here, we report the near-complete backbone and partial side-chain NMR resonance assignments of the C-terminal domain of human 5MP1 (residues 250-419), carrying a W404E substitution that disrupts dimerization. The WT protein forms a dimer at NMR concentrations, which increases the effective size of the protein and also causes disappearance of peaks corresponding to aminoacids at the dimer interface due to conformational exchange. Backbone resonance assignments were completed for 96.4% of the non-proline residues. Secondary structure was analyzed using Chemical Shift Index (CSI) and compared with the AlphaFold structural model. Regions of disagreement between the experimental and computational secondary structure assignments were further examined using 15N-NOESY-HSQC spectra, allowing experimental validation of local structural features. While the AlphaFold model accurately reproduces the overall fold of the 5MP1 C-terminal domain, several localized discrepancies were identified, particularly near the N- and C-terminal regions of the domain, where experimental NMR data support alternative secondary structure assignments. These resonance assignments and experimentally validated structural features provide a foundation for future investigations of the molecular interactions, dynamics, and functions of 5MP1 in translation initiation.

Cisplatin resistance remains a major barrier in head and neck squamous cell carcinoma (HNSCC) treatment. ATP-binding cassette (ABC) transporters contribute to chemoresistance by limiting intracellular drug accumulation. Bitter taste receptor 10 (T2R10) has been implicated in ABC transporter regulation, but its role in HNSCC remains undefined. HNSCC cell lines were treated with T2R10-agonist caffeine (100 or 200 M), cisplatin, or a combination, and viability was assessed by crystal violet assay. T2R10 promoter activity and expression following caffeine exposure were evaluated using a promoter-driven mCherry reporter and RT-qPCR. ABC transporter expression was measured after caffeine treatment and T2R10 gene (TAS2R10) knockdown or overexpression. Associations between tumor TAS2R10 expression and survival were assessed using TCGA data through GEPIA2. Caffeine enhanced the cisplatin-associated reduction in viability in a cell line- and concentration-dependent manner, with the strongest effect seen in UM-SCC47. A significant effect was observed in FaDu at 200 M of caffeine, and minimal response in RPMI 2650. RPMI 2650 cells and FaDu cells exhibited lower baseline TAS2R10 expression and RPMI 2650 cells did not demonstrate enhanced cisplatin sensitivity following caffeine treatment. Caffeine treatment increased TAS2R10 promoter activity and expression and was associated with decreased ABCG2 expression. TAS2R10 knockdown increased ABCG2 and ABCF1 expression, whereas TAS2R10 overexpression reduced ABCG2 and ABCC1 expression. High tumor TAS2R10 expression was associated with improved disease-free survival (log-rank p=0.0071; HR=0.61) but not overall survival. Caffeine enhances cisplatin sensitivity in selected HNSCC models. Caffeine exposure is associated with increased TAS2R10 expression and reduced expression of chemoresistance-associated transporters, particularly ABCG2.

Hurthle cell carcinoma (HCC) is an aggressive form of thyroid cancer. While mitochondrial DNA mutations and chromosomal losses have been identified in HCC, isoform switching, and its functional consequences remain uncharacterized. This study reanalyzed NCBI GEO dataset GSE228870 (n = 32), using Salmon and IsoformSwitchAnalyzeR() to identify isoform switching. The analysis resulted in 371 switches across 335 genes showing functional consequences including loss of protein domains, shorter open reading frames (ORFs), loss of signal peptides and novel sub-cellular localizations. Most significant isoform switches (q-value < 0.05, |dIF| > 0.1) were observed in LAMA2, LSP1, MAD2L2, FBLN2 and CXCL12, implicating extracellular matrix dysregulation, DNA damage response, immune signaling and cytoskeleton regulation. These genes are expressed in normal thyroid (median TPM 20.69, 11.66, 14.79, 134.1 & 80.76). However, specific isoforms of LAMA2 and MAD2L2 are not expressed in normal thyroid, explaining tumor-specific expression in HCC. Alternative transcription termination site (ATTS) gain was significant, suggesting altered 3 end in HCC transcripts. TCGA SpliceSeq showed LSP1, FBLN2 and CXCL12 undergo alternative promoter (LSP1 exon1 PSI=94.5%, FBLN2 exon2 PSI=99.0%) and alternative termination (CXCL12 exon3.3 PSI=53.9%) in thyroid cancer, suggesting ATTS and alternative transcription start site (ATSS) as shared splicing dysregulation mechanisms. This is the first systematic characterization of isoform-level dysregulation in HCC.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

In TMPRSS2 fusion-positive prostate cancer, ERR is involved in regulation of ERG and promotes the androgen receptor independent signaling in the cancer progression. The ERR binds to the ERREs (estrogen-related receptor response elements) present at -5042 bp of the TMPRSS2- promoter and enhances the ERG overexpression that causes prostate cancer progression. To dissect the structural basis of the ERR recognition to the TMPRSS2 promoter DNA, we have purified the full-length ERR (ERRFL), NTD deleted construct (ERR{Delta}NTD), and the DNA-binding domain (ERRDBD) proteins and performed the binding analysis with 30 bp TMPRSS2-promoter DNA (5' -AGTCCAAGGTCGGTGGATC ACAAGGTCAGG-3'). Circular dichroism analysis showed that all three ERR proteins adopt native secondary structures. DNA binding induced subtle changes in the secondary structures, while enhancing the thermal stability (Tm) of all ERRa proteins. Binding analysis showed that ERRDBD bound weakly to the DNA, whereas ERRFL and ERR{Delta}NTD exhibited substantially higher affinities ~120-fold and ~131-fold than ERRaDBD, respectively. Small-angle X-ray scattering (SAXS) analyses revealed a dimeric ERRFL structure and an ERRFL-DNA complex (2:1) structure in solution and fitted well with Alpha Fold model of apo and DNA bound complex of ERRFL. Furthermore, 100 ns dynamics simulations on apo and DNA-bound ERRa proteins showed that all proteins remained structurally stable, with flexibility largely confined to loop regions of ERRa proteins. Our biophysical, DNA binding and structural analyses have revealed the mechanism involved in ERR recognition of the TMPRSS2- promoter DNA, which provides insight into ERR-mediated transcriptional regulation and development of anticancer drugs against ERR-driven prostate cancer.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

Oncogenic mutations are known to occur at non-random rates and specific hotspots, yet the structural and thermodynamic factors underlying these hotspots remain poorly understood. This study presents molecular dynamics simulations of 22 cancer hotspots across 10 oncogene families, simulated as histone-DNA complexes. Across 18 of 22 mutations, thermodynamic effects were predominantly localized to within 6 [A] of the DNA-histone interface (mean capture 104%), with van der Waals interactions driving the effect at the thermodynamic extremes, indicating that oncogenic mutations alter nucleosome binding through precise local contact changes rather than global structural rearrangements. Further, analysis revealed a gene-family correlated pattern in thermodynamic stability. RAS family mutations showed a consistent trend toward nucleosome stabilization (mean {Delta}{Delta}G = -4.82 kcal mol-1), while kinase domain mutations trended toward destabilization (mean {Delta}{Delta}G = +40.12 kcal mol-1), a difference reaching statistical significance in this exploratory analysis (Mann-Whitney U, p = 0.008). This interface-specific mechanism, combined with the gene-family-correlated thermodynamic pattern, provides a biophysical framework for understanding nucleosome-level contributions to cancer hotspot mutation biology.

Three-dimensional (3D) tumor models exhibit drug responses that differ from conventional 2D cultures. However, how cellular metabolism dynamically evolves across culture models and during drug treatment remains poorly understood. We compared the responses of MCF-7 cells in 2D and 3D environments to 5-fluorouracil (5-FU), combining {superscript 1}H NMR cellular metabolomics with viable cell counts, the GLUT1-positive population, gene expression, ATP activity, and morphology. Principal component analysis revealed that culture dimensionality, rather than 5-FU treatment, was the primary driver of metabolic flux variation. 3D spheroids exhibited higher glycolytic flux at 72h. Importantly, this elevated glycolysis reflected a higher per-cell flux in larger spheroids rather than an increased cell number. We further observed a higher proportion of GLUT1-positive cells and increased HK2 expression in 3D culture, together with an epithelial phenotype characterized by increased CDH1 and decreased VIM expression. Functionally, 3D displayed maintaining higher cell viability, ATP activity following treatment. Together, these findings suggest that 3D architecture promotes a metabolically defensive phenotype and cellular metabolic behaviors is associated with morphology, which may inform future drug screening model selection. BlurbTime-resolved NMR metabolomics reveals that 3D culture architecture, rather than 5-FU treatment, defines the metabolic phenotype of MCF-7 cells, linking spheroid morphology with per-cell glycolytic activity, ATP preservation, and reduced chemotherapy sensitivity. Synopsis O_FIG O_LINKSMALLFIG WIDTH=197 HEIGHT=200 SRC="FIGDIR/small/740505v1_ufig1.gif" ALT="Figure 1"> View larger version (41K): org.highwire.dtl.DTLVardef@cf6b74org.highwire.dtl.DTLVardef@197a4aeorg.highwire.dtl.DTLVardef@bcb5f8org.highwire.dtl.DTLVardef@14ec1ac_HPS_FORMAT_FIGEXP M_FIG C_FIG Graphical abstract -- Time-resolved metabolic and phenotypic responses of 2D and 3D MCF-7 cultures to 5-fluorouracil. Bullet pointsO_LITime-resolved {superscript 1}H NMR cellular metabolomics combined with multilevel phenotypic readouts reveals that culture dimensionality, rather than 5-FU chemotherapy, is the dominant determinant of metabolic phenotype in MCF-7 breast cancer cells. C_LIO_LI3D spheroids develop a glycolysis-dominant metabolic state with markedly elevated GLUT1 populations and HK2 transcript, accompanied by substantially reduced sensitivity to 5-FU compared with 2D monolayers. C_LIO_LISpheroid morphology correlates with cellular metabolic flux: larger and more elongated spheroids exhibit higher glycolytic activity on a per cell, independent of spheroid cell number. C_LIO_LIThe 3D-defined metabolic state buffers ATP under 5-FU stress and reinforces an epithelial transcriptional program (CDH1{uparrow}, VIM{downarrow}), arguing that culture architecture must be considered when interpreting preclinical drug responses. C_LI

In terms of cancer-related death, non-small cell lung cancer (NSCLC), the worlds leading cause, highlights the need for continued research into the genetic factors that influence tumor growth. Long non-coding RNAs (lncRNAs) are now well recognized as essential regulators of oncogenic signaling cascades; nevertheless, the specific role and molecular basis of the SOX2 overlapping transcript (SOX2OT) in NSCLC are not entirely understood. This study examined the functional importance of SOX2OT and its regulatory interactions with tumor-suppressive microRNAs in NSCLC cells. In A549 and Calu-3 cells, RNA interference-mediated SOX2OT silencing dramatically reduced cellular proliferation, migration, and invasiveness. Moreover, SOX2OT knockdown was associated with inhibition of epithelial-mesenchymal transition (EMT), alongside induction of cell cycle arrest and activation of apoptotic pathways. Integrated transcriptomic profiling and bioinformatic prediction analyses identified miR-143 as a putative downstream effector of SOX2OT activity. Consistently, depletion of SOX2OT resulted in marked elevation of miR-143 expression, which corresponded with downregulation of oncogenic mediators, including STAT3, EZH2, and CXCL13. As a result of SOX2OT suppression, both the transcript and the protein levels of PTEN were restored. Further functional characterization demonstrated that SOX2OT knockdown inhibits EMT progression by decreasing mesenchymal markers and EMT-related transcription factors (TFs) while concomitantly enhancing epithelial marker expression. Collectively, these findings suggest that SOX2OT contributes to NSCLC pathogenesis through regulation of a miR-143-centered signaling network that influences oncogenic signaling, cellular survival, and metastatic potential. Targeting the SOX2OT/miR-143 regulatory axis may therefore represent a promising therapeutic approach for NSCLC, while also underscoring the broader importance of lncRNA-mediated post-transcriptional regulation in lung cancer biology.

The TMEM16 family of calcium-activated phospholipid scramblases (CaPLSs) and chloride channels (CaCCs) performs diverse physiological functions that include regulation of blood coagulation and apoptotic signaling, through a shared ten-transmembrane-helix (TM) architecture organized around a hydrophilic lipid-translocating groove. Mechanistic studies of TMEM16 family members have been hampered by the absence of a unified positional reference framework that would permit direct comparison of structurally equivalent residues across paralogs with different sequence numbering systems. Here we introduce a generic numbering scheme for TMEM16 scramblases (GNS-TMEM16), modeled on the Ballesteros & Weinstein system established for class A G protein-coupled receptors. A reference alignment (TMEM16-RA) was constructed from twelve human and mouse TMEM16 scramblases (TMEM16C/D/E/F/G/J) using structure-based ClustalW alignment of the ten TM helices. From this alignment, a TM-specific reference residue (TsRR) was identified for each helix by hierarchical application of three criteria: (1) 100% conservation in the core TMEM16-RA; (2) conservation in an augmented reference alignment (TMEM16-ARA) incorporating a group of phylogenetically more distant homologs composed of nhTMEM16, afTMEM16, TMEM16K, TMEM16A, and TMEM16B; and (3) structural and functional considerations, including helix-perturbing character, groove localization, conserved motif membership, and central TM position. The resulting ten TsRRs are Y1.50, W2.50, R3.50, E4.50, F5.50, P6.50, E7.50, D8.50, W9.50, and E10.50, and are illustrated in mTMEM16F. Each residue is assigned the identifier N.m(k), where N is the TM number, m is the position relative to the TsRR (for which m = 50), and k is the absolute sequence number. Loop residues receive dual identifiers referenced to the TsRRs of both flanking helices. Application of the GNS-TMEM16 is illustrated with the comparisons of the groove-opening measurements using pairwise distances between residues identified by their N.m indices to be corresponding across mTMEM16F, afTMEM16, and nhTMEM16. The results bring to light the advantages of corresponding residues identification in different TMEM16 proteins and show that the mammalian scramblase undergoes substantially larger separation at the extracellular groove entrance than either fungal homolog. Comparison of mutagenesis data guided by N.m correspondence shows at the conserved (E3.55,R6.26) salt-bridge locus, Ala substitution reduces activity more than 100-fold in nhTMEM16 but less than 2-fold in afTMEM16, illustrating that the GNS identifies structural equivalence of position without implying functional equivalence of the residue, which is a distinct advantage of GNS in providing mechanistic interpretation across paralogs. Also described is a protocol for extending the GNS-TMEM16 to uncharacterized protein sequences, including AlphaFold-predicted models, using structural superposition to mTMEM16F. Thus, the presented GNS-TMEM16 provides a stable positional reference for the integration and comparative analysis of structural, computational, and functional data across the TMEM16 family, utilizing a construction strategy applicable to yet other polytopic membrane protein families sharing a common transmembrane fold.

Bis(monoacylglycero)phosphate (BMP) is a signature lysosomal phospholipid that supports the catabolic functions of the lysosome. We recently demonstrated that BMP deficiency is associated with a variant of Batten disease, a neurodegenerative lysosomal storage disorder. This observation led us to investigate how BMP deficiency contributes to the previously reported accumulation of the ATP synthase c-ring in the lysosomes of Batten disease patients and preclinical models. The c-ring is an inner mitochondrial membrane (IMM) protein complex that interacts with cardiolipin, a mitochondrial phospholipid that shares structural features with BMP. Based on this, we hypothesised that BMP may perform an analogous function to cardiolipin in lysosomes. Specifically, we proposed that BMP preferentially interacts with the c-ring, dispersing it within lysosomal membranes and facilitating its degradation. To test this hypothesis, we conducted all-atom molecular dynamics simulations to examine the interactions of various BMP variants with the c-ring of human ATP synthase under different membrane conditions. We observed leaflet-specific preferential interactions of BMP with the protein interface. Replacement of BMP with anionic POPG lipids resulted in a lower binding affinity to the c-ring, indicating the importance of BMPs unique structure with respect to protein binding. Furthermore, BMP enrichment was enhanced when using the physiologically relevant di-22:6 BMP variant in membranes containing polyunsaturated lipids and cholesterol. Overall, our study suggests that BMP promotes lysosomal c-ring degradation via c-ring co-localisation, whereas BMP deficiency in Batten disease drives c-ring accumulation. Significance statementBatten disease comprises a group of neurodegenerative lysosomal storage disorders that primarily affects infants and children. Several recently identified disease genes are linked to the biosynthesis of lipid bis(monoacylglycero)phosphate (BMP), a lysosomal lipid required for normal lysosome function. However, how BMP deficiency contributes to disease pathology remains poorly understood. Here, we investigate the relationship between BMP deficiency and lysosomal accumulation of the ATP synthase c-ring, a hallmark pathological feature of Batten disease. Using molecular simulations, we show that BMP preferentially interacts with the c-ring, suggesting that BMP deficiency promotes c-ring aggregation within lysosomes. Our study provides a mechanistic starting point for understanding how BMP deficiency disrupts lysosomal function in Batten disease

Mutations in mitochondrial DNA (mtDNA) compromise ETC activity and impair oxidative phosphorylation. Since eukaryotic cells contain multiple copies of mtDNA, the resulting phenotype depends on the proportion of mutant mitochondrial genomes (the heteroplasmy level). Using isogenic cell lines carrying similar mtDNA deletions, a linear decline in cellular respiration was observed as mitochondrial DNA heteroplasmy increased. Despite this, cellular redox imbalance did not change until heteroplasmy exceeded 50%. As heteroplasmy increased past 70%, cells also exhibited an integrated stress response (ISR) and impaired translation was observed. These defects were reversed by either addition of asparagine or overexpression of pyruvate carboxylase (PC). The dependence on exogenous asparagine in other respiration-deficient cells was found to correlate inversely with the PC expression level. For example, patient-derived thyroid tumor cells, harboring high heteroplasmy for a Complex I mtDNA mutation and low levels of PC, exhibited asparagine auxotrophy, and L-asparaginase treatment suppressed tumor growth. Together, these findings demonstrate a role for mitochondrial pyruvate carboxylase in cellular asparagine synthesis under conditions of compromised respiratory activity.

P-type ATPases maintain cellular metal homeostasis by coupling ATP hydrolysis to the active translocation of ions across biological membranes. Within this superfamily, P1B -type ATPases play a central role in the transport and detoxification of transition metals such as Cu+ and Ag+. Despite the availability of several crystallographic structures for bacterial ATPases, molecular determinants governing ion selectivity, coordination dynamics, and transport efficiency are not yet fully resolved, particularly for plant homologs, whose conformational landscapes have not been explored in detail at atomistic resolution. In this study, we combine molecular dynamics (MD) simulations with enhanced sampling through metadynamics to characterize the conformational dynamics and ion-binding mechanisms of representative heavy-metal P-type ATPases. By comparing bacterial and plant homologs in phosphorylated and dephosphorylated states, we examine how phosphorylation is associated with changes in structural plasticity, channel accessibility and ion coordination environments. The simulations indicate phosphorylation-dependent rearrangements of the cytosolic domains and transmembrane helices, which may alter access to the metal-binding pathway. Enhanced sampling identifies distinct conformational basins and possible transition pathways along functionally relevant collective variables, allowing characterization of the free-energy landscape associated with ion displacement within the transmembrane region. Together, these results provide an atomistic perspective on the relationship between metal coordination and protein conformational dynamics. Overall, our findings support a model in which phosphorylation state, structural dynamics, and local coordination environments collectively contribute to the differential behavior of metal ions in P1B -type ATPases.

The androgen receptor is a hormone-dependent transcription factor that regulates a wide range of physiological processes and plays a pivotal role in the development of prostate cancer. Its 555-residue, intrinsically disordered, N-terminal domain is involved in the modulation of transcriptional activity by recruiting co-regulators and mediating the formation of biomolecular condensates. This study reports on the characterisation of a conserved domain located in the C-terminal region of the androgen receptor N-terminal domain, at atomic level, using nuclear magnetic resonance spectroscopy. This proline rich region exhibits extensive conformational heterogeneity driven by highly populated cis proline conformers that are stabilised through interactions with adjacent aromatic residues. We demonstrate that the cis-proline population is modulated by phosphorylation as well as cancer-associated mutations. This suggests that proline driven conformational heterogeneity at the C-terminal region of androgen receptor N-terminal domain is involved in the regulatory function of this transcription factor.